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Abstract 


This  Report  discusses  the  special  relativistic  and  general  relativistic 
effects  on  electromagnetic  propogation  for  artificial  satellite  applications. 
In  particular  the  rigorous  special  relativistic  Doppler  shift  between  the 
received  and  transmitted  frequencies  is  derived.  The  two  radars  need  not  be 
colocated  and  are  assumed  not  to  be  in  the  same  inertial  reference  frame. 
Numerical  estimates  of  various  approximations,  second  order  series  approxima¬ 
tions,  and  lowest  order  post-Newtonian  general  relativistic  terms  are  all 

performed.  This  work  was  spurred  by  the  observational  determination  of  the 
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I.  BASIC  FORMULAS 

The  central  formula  from  the  theory  of  special  relativity  pertinent  to 
this  discussion  is  the  Doppler  shift  equation.  This  result  may  be  derived  in 
a  variety  of  ways  and  may  be  looked  upon  as  an  expression  of  the  coordinate 
system  invariance  of  the  phase  of  an  electromagnetic  wavetrain.  As  the  phase 
is  dependent  upon  the  scalar  product  between  two  4- vectors,  we  are  used  to 
seeing  the  Doppler  shift  results  displayed  as  two  equations— one  from  the 
temporal  part  of  the  4- vector  scalar  produce  and  one  from  the  spatial  part. 
Actually,  because  of  the  nature  of  the  Lorentz  transformation  for  a  4-vector, 
the  Doppler  shift  formulas  can  be  broken  into  three  parts. 

Consider  two  inertial  reference  frames  in  relative  motion.  Let  their 
relative  velocity  be  given  by  the  3-vector  v^  =  cB^  where  c  is  the  speed  of  light. 
Define  y  as 

y  -  d-e2)  -1/2 

Then  the  Lorentz  transformation  of  an  arbitrary  4-vector  A  =  (Aq,A-j  .Ag.A.^)  = 
(Aq,A)  is  given  by 

a0  =  y^A0  " 
a„  =  y(A(|  -  SA0) 


where  one  observer  measures  A  =  (A^.A)  and  the  other  inertial  observer  per¬ 
ceives  a  *  ( »a.) •  The  notation  ii  and  l  refer  to  directions  parallel  to  and 
perpendicular  to  v_. 
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For  a  plane  wave  of  frequency  F  and  wave  vector  J<  the  result  Is 
f  =  y(F  -  i-K) 
kn  =  y(K„  "  &F) 


Because  of  the  special  nature  of  light  waves  in  special  relativity  these  may 
be  more  simply  written  as  (F=Kc) 

f  =  yF(1  -  Bcose) 


tan  6 


sin  0 
y(cos  Q-e) 


(!) 


where  G  (0)  is  the  angle  between  J<  (k.)  and  v.  The  inverse  to  Eqs.  (1)  may  be 
obtained  by  interchanging  lower  case  and  upper  case  letters  and  changing  the 
sign  of 

Another  way  to  describe  the  meaning  of  the  angle  0  (0)  is  that  if  the 
inertial  observer  (the  inertial  observer  moving  with  velocity  v)  perceives 
the  direction  of  propogation  of  the  ray  to  be  along  the  unit  vector  (nj, 
then  6  cos  O  =  J3*N^  (Scos  0  =  6*nJ.  In  this  notation 


f  =  yF(1  -  M) 


■  (F/y)/(1  +  £-n) 


(2) 


The  second  of  Eqs.  (1)  is  an  aberration  effect  that  is  present  even  if  the  two 
observers  are  instantaneously  colocated. 


II.  THE  SCENARIO 

The  transmitting  radar  emits  a  pulse  of  radiation  out  into  space.  This 
transmitted  pulse  (T  subscript  below)  is  reflected  by  an  artificial  satellite 
sometime  later.  Even  later  the  reflected  pulse  is  intercepted  by  the  receiv¬ 
ing  radar.  See  Fig.  1.  The  receiving  radar  may  be  colocated  with  the  trans¬ 
mitting  radar  (possibly  identical  with  it)  or  not.  The  transmitting  radar 
knows  the  frequency  of  the  transmitted  pulse  (fy)  as  measured  in  the  inertial 
reference  frame  which  is  instantaneously  comoving  with  it  at  the  time  of 
transmission.  Continuing  in  this  anthropomorphic  vein,  the  transmitting  radar 
also  knows  the  direction  of  propogation  of  the  outgoing  pulse  as  perceived  in 
this  same  inertial  frame  (unit  vector  =  rty). 

The  receiving  radar  has  available  to  it  similar  knowledge  concerning  the 
received  pulse  (R  subscript  below).  In  particular  this  radar  observes  an 
incoming  pulse  of  frequency  fR  coming  from  a  direction  n^.  Both  radars  are 
also  aware  of  their  velocities  (at  the  time  of  transmission  and  reception 
respectively)  relative  to  an  inertial  observer  located  at  the  center  of  the 
Earth.  Indeed,  the  geocentric  inertial  observer  perceived  the  transmitting 
radar,  instantaneously  moving  with  a  velocity  Vy  =  cj3y,  to  transmit  a  pulse 
of  frequency  Fy  in  a  direction  Ny.  Sometime  later  an  artificial  satellite 
with  velocity  =  c6g  (relative  to  the  geocentric  inertial  frame)  reflected 
this  pulse.  Later  still  the  geocentric  inertial  observer  perceived  the 
intercepting  radar,  instantaneously  moving  with  a  velocity  v^  =  cg^,  to 
receive  a  pulse  of  frequency  FR  from  the  direction  N^. 

The  desideratum  in  this  problem  is  to  express  the  frequency  component 
of  the  Doppler  shift,  Af  =  fR-fy,  in  terms  of  known  (the  velocities)  or 
observed  (the  directions  and  the  frequencies)  quantities. 
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Fig.  1.  Schematic  representation  of  the  transmitted  (T)  and  received  (R) 
radar  signals.  The  displacement  may  be  due  to  the  Earth's  rotation  or 
the  non-colocated  radars. 


III.  RIGOROUS  SPECIAL  RELATIVISTIC  RESULT 


Define  yR,  yt  by 


Y,  ■  0  -  «S> 

YT  *  0  -  6?)  -1/2 


Then  from  Ea.  (2)  the  relationship  between  the  geocentric  inertial  observer's 
perception  of  the  frequency  of  the  transmitted  pulse  (Fj)  and  the  frequency 
detected  by  an  observer  instantaneously  comoving  with  the  transmitting  radar 
(fT)  is 


fT  =  Fy  Yy  0  -  Sy-Ny) 

=  (Fy/Yy)/(1  +  By’Hy) 


(3) 


The  next  step  in  the  scenario  is  the  reflection  of  this  outgoing  pulse  by 
an  artificial  satellite.  In  the  geocentric  inertial  reference  frame  the  arti¬ 
ficial  satellite  has  a  velocity  =  cBg  and  reflected  the  transmitted  pulse 
in  the  direction  N^.  An  observer  instantaneously  comoving  with  the  artificial 
satellite  perceived  the  transmitted  pulse  approaching  with  frequency  fs> 

fS  =  FT  yS  ^  "  — S*— T^ 

where  ys  =  (1  -  3^)  ~^2.  The  reflected  pulse  receded  from  the  satellite's 
surface,  relative  to  the  geocentric  inertial  observer,  at  a  frequency  FR. 


This  frequency  Is  related  to  f^  via 
fs  =  Fr  YS  (1  -  85*!^) 

So,  utilizing  the  last  two  formulas,  I  can  write 
Fr  ■  (fs/Ys)/(l  -  Is’Nr) 

-  Ft  (1  -  4-Nt)/(1  -  ij-N,,)  (4) 

It  is  important  to  note  that  the  relationship  between  the  received  and  trans¬ 
mitted  pulse  frequencies  (in  the  geocentric  inertial  frame)  contains  two  extra 
special  relativistic  components.  One  is  equal  to  ys»  the  other  to  1/ys. 

Hence  their  cumulative  effect  is  nil  and  the  final  relationship  has  the  appear¬ 
ance  of  a  Newtonian  formula. 

Finally  the  receiving  radar  detects  a  pulse  of  frequency  fR  coming  from 
the  direction  at  the  instant  when  its  velocity  is  v^  =  cj^  relative  to  the 
geocentric  inertial  observer.  By  using  Eqs.  (2),  (4),  and  then  (3)  fR  can  be 
expressed  as 

fR  =  FR  yR  ^  ~  V 

=  (VtrJ/O  + 

=  Ft  Yr  (1  -  6^)0  -  6s*Nt)/(1  -  6S%)  (5) 


=  fT  (yo/Yt)0  -  1d*Nd)(1  ‘  Bc-Nt) 


Lastly  the  Doppler  shift  in  frequency,  Af  =  fR  -  fy,  can  be  expressed  as 


Af / f T  = 


i  -4Y/2  /i 


i  -  e; 


1  -  it-*t 


1  £R 
1  '  — S'^R 


Note  that  Af  is  the  difference  between  two  frequencies  measured  in  two  differ¬ 
ent  inertial  frames  and  that  Af/fy  in  Eq.  (6)  is  a  quantity  expressed  in  terms 
of  velocities  and  directions  relative  to  the  geocentric  inertial  reference 
frame. 

* 

Irwin  Shapiro  discussed  a  similar  problem  in  the  interplanetary  context. 
In  his  case  the  "artificial  satellite"  was  Venus,  the  transmitting  and  receiv¬ 
ing  "radars"  were  the  Earth  at  the  instants  of  transmission  and  reception, 
and  the  ultimate  reference  frame  was  a  heliocentric  one  instead  of  a  geocentric 
one.  To  clearly  see  the  correspondence  between  Shapiro's  result  and  Eq.  (6), 
let  T  1  ,\r  -+3,  S  -+2,  Ny  ey2>  and  -*  e^. 

Finally,  when  the  receiving  radar  and  the  transmitting  radar  are 
colocated,  BR  =  By  (but  3^  t  By).  Hence  Eq.  (6)  reduces  to 

Af/fy  =/- - &AV- - -  1  (7) 

V  “  — T*— T/  Y  “  h'\) 

Once  again,  while  this  formula  has  the  appearance  of  one  which  conforms  to 
Galilean  relativity,  it  was  derived  (in  full  generality)  within  the  context 


of  a  more  correct  physical  theory  of  nature.  Its  final  simplicity  is  due 


I.  I.  Shapiro,  Bull,  Astr.  XIV,  201,  1965. 


to  a  degenerate  geometry  (transmitting  radar  s  receiving  radar)  and  the 
fortuitous  cancellation  (by  division)  of  two  of  the  four  special  relativistic 
terms  ■‘mpllclt  In  Eq.  (7). 
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IV.  APPROXIMATIONS 


Since  the  artificial  satellite  is  nearby  (or  light  travels  fast)  an 
approximation  that  appears  tempting,  when  the  transmitting  and  receiving 
radars  are  identical,  is  to  ignore  the  displacement  of  the  radar  during  the 
pulse  travel  time  (radar-satellite-radar).  The  round  trip  time  is  =*0^3  even 
for  a  near-stationary  satellite  and  the  Earth's  rotation  rate  is  slow  (=*15"/ 
sec  so  neglecting  the  displacement  is  an  error  of  2  parts  in  105}.  If  the 
radar  doesn't  move  (or  the  pulse  is  infinitely  fast)  then  it  follows  that 
the  directions  of  transmission  and  reception  are  antiparallel.  Hence, 

N-j.  -  -N^  =  N  and  =  -n^  =  n.  Furthermore,  if  the  radar's  displacement  is 
being  neglected,  then  the  approximation  that  jiy  =  =  8  (say  at  the  instant 

of  reflection  by  the  satellite)  can  also  be  justified.  Now 


Af/fT 


1  +  I*N  \ 

}  *  HV  " 1 


(8) 


Note  that  ignoring  the  displacement  of  the  radar  during  the  time  of 
flight  of  the  pulse  is  not  the  same  as  ignoring  the  relative  motion  of  the 
radar  with  respect  to  the  geocentric  inertial  observer.  Were  I  to  do  that 
then  I  could've  set  £  =  0  in  Eq.  (8).  This  is  incorrect  and  would  lead  to 
the  neglect  of  an  aberration  effect  [see  Eq.  (1)].  This  aberration  effect 
is  distinct  from  any  parallactic  displacement  and  solely  results  from  the 
relative  motion  of  the  two  observers.  Moreover,  as  a  glance  at  Eq.  (1)  will 
show,  it  is  a  first  order  effect  (e.g.,  it  depends  on  6). 

Since  8  and  Bs  are  small  (<10~5)  a  power  series  expansion  of  Eq.  (8) 
would  seem  to  be  appropriate.  The  two  leading  terms  are 


9 


(9) 


Af/fT  a  -2  A$*N  +  2  (A_B*N)2 

where  A3  =  B.s  -  8..  AMs  the  relative  velocity  of  the  artificial  satellite 
with  respect  to  the  radar  as  measured  by  the  geocentric  inertial  observer. 

A£*N^  is  the  radial  velocity  of  the  satellite,  relative  to  the  radar,  as  per¬ 
ceived  by  the  geocentric  inertial  observer.  This  is  because  AjS  is  not  the 
relative  satellite-radar  velocity  in  the  radar's  comoving  reference  frame  nor 
is  the  direction  to  (the  "line-of-sight"  of)  the  satellite  in  the  radar's 
comoving  reference  frame  (]i  is).  The  correct  radial  velocity  of  the  satellite 
relative  to  the  radar,  in  the  radar's  comoving  reference  frame,  can  be 
obtained  by  performing  a  Lorentz  transformation  of  A6.  and  then  forming  the 
scalar  product  of  the  transformed  vector  with  n.  AJ3  is  the  Galilean  result 
for  the  relative  velocity  of  the  satellite  with  respect  to  the  radar.  Finally 
note  that  while  the  difference  between  Aj3  and  the  relative  velocity  of  the 
artificial  satellite  with  respect  to  the  radar  (in  the  radar's  comoving 
reference  frame)  is  of  second  order,  the  difference  between  and  n^  is  first 
order.  As  the  point  of  this  derivation  is  the  full  set  of  second  order  terms 
neither  of  these  effects  can  be  neglected.  [If  you  object  that  only  appears 
in  Eq.  (9)  dotted  into  AB-,  and  that  A3  is  of  first  order,  the  net  effect  of 
N  f  in  is  of  the  second  order  In  Af/fy.  The  aberration  effect  is  still  of  the 
first  order.] 

Equation  (9)  Is  not  in  a  form  suitable  for  computation  because  N_  is  not 
observed.  The  relative  velocity  A3.  *  £s  -  £  Is  expressed  geocentrically  and 
this  is  the  coordinate  system  wherein  It  is  simplest  to  compute  both  £  (the 
radar's  geocentric  velocity)  and  3g  (the  artificial  satellite's  geocentric 
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velocity).  If  I  return  to  Eq.  (6)  and  use  Eq.  (2),  then  two  of  the  geocentric 
directions  can  be  eliminated  yielding 


Af/fT  =  (yt/yr) 


-  1 


(10) 


The  only  way  I  can  see  to  eliminate  the  other  two  geocentric  inertial  unit 
vectors  is  to  replace  them  with  unit  vectors  in  the  artificial  satellite's 
comoving  reference  frame.  As  this  doesn't  facilitate  matters,  I  have  demurred. 

It  is  still  true  that  By  =  8R  if  the  receiving  radar  is  colocated  with 
the  transmitting  radar,  whence  Yy  =  YR*  If  in  addition  I  ignore  the  displace¬ 
ment  of  the  radar  during  the  pulse's  transit  time,  then  Ny  =  =  N  and 

rvj.  =  -n^  =  n^  again.  Af/fy  reduces  to 


=  -2  Bs-N  +  2  B*n  -  4  (Bs'N)(B*n)  +  2  (6-n)2 
+  2  (Bs-N)2 


Taking  both  B  and  8S  to  be  comparably  small  (in  practice  8$  _>  3  perhaps 
10B)  I  can  replace  N  by  n.  in  any  quadratic  term.  After  making  this  approxi¬ 
mation  I'm  left  with  Eq.  (12), 

Af/fy  =  -  2  BS*N  +  2  B*n  +  2  (A3*n)2  (12) 
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Finally,  I'm  going  to  write  £  as  n.  +  6n_  where  the  norm  of  6n  is  of  order  6. 

Then  the  expression  for  Af/fy  reduces  to  (A3.  =  £s  -  S  still) 

Af/fy  »  -  2  AS-n  +  2  (A£*n)2  -  2  3$*  6n  (13) 

The  difference  between  the  expression  in  Eq.  (13)  and  in  Eq.  (9)  clearly  shows 
the  effect  of  aberration. 

A0*n  is  still  not  the  radial  velocity  of  the  artificial  satellite  with 
respect  to  the  radar  (in  any  reference  frame)  because  A£  is  still  the  relative 
velocity  of  the  satellite  with  respect  to  the  radar  In  the  geocentric  Inertial 
frame.  The  fact  that  n^  is  the  direction  to  the  artificial  satellite  in  the 
radar's  comoving  inertial  frame  makes  A6*n  a  better  approximation  to  the  rela¬ 
tive  radial  velocity  than  Is  Ajj-f^.  To  complete  the  computation  a  Lorentz 
transformation  of  AB  would  be  required. 

If  Eq.  (13)  or  Eq.  (9)  Is  viewed  as  the  ultimate  result  of  the  analysis 
then  the  reader  should  be  explicitly  aware  of  two  implicit  assumptions.  One 
concerns  the  special  theory  of  relativity,  the  other  the  general  theory.  I 
have  spoken  above  of  a  geocentric  Inertial  observer.  The  implication  of  such 
a  phrase  is  that  relative  to  any  Inertial  reference  frame  the  velocity  of  the 
geocenter  has  been  constant  throughout  the  pulse  travel  time  (radar-satellite- 
radar).  As  the  Earth  revolves  about  the  Earth-Moon  barycenter,  and  the  Earth- 
Moon  barycenter  revolves  about  the  solar  system's  barycenter,  and  the  Sun  is 
accelerated  relative  to  the  local  standard  of  rest,  and  the  local  standard  rest 
revolves  about  the  galactic  center,  .  .  .  this  hypothesis  is  known  to  be  false. 
The  numerical  consequences  of  such  an  assumption  will  be  treated  in  the  next 
Section. 
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The  general  relativistic  effect  not  obviously  Included  concerns  the 
change,  if  any,  between  the  gravitational  potentials  of  the  transmitting  and 
receiving  radars.  As  the  practical  case  has  these  colocated,  the  essence  of 
the  approximation  Is  to  ignore  the  motions  of  the  other  bodies  in  the  solar 
system  during  the  round  trip  flight  time  of  the  pulse.  This  too  needs  to  be 
numerically  Investigated. 
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V.  NUMERICAL  CONSIDERATIONS 

There  are  actually  two  general  relativistic  effects  that  can  yield 
systematic  sources  of  error  in  the  Newtonian  reduction  of  radar  observations. 
The  first  is  the  general  relativistic  Doppler  shift  mentioned  above.  The 
amplitude  of  the  change  is  proportional  to  the  difference  between  the  gravi¬ 
tational  potentials  of  the  transmitting  and  receiving  radars  (in  the  lowest 
order  post-Newtonian  approximation).  In  particular  the  shift  in  wavelength 
AX  corresponding  to  the  frequency  shift  Af  is  given  by 

AX/X  =  A4/c2 

where  is  the  negative  of  the  usual  Newtonian  gravitational  potential.  For 
a  particle  on  the  Earth's  surface  *  GM^/R^  but  this  expression  ignores  both 
the  oblateness  of  the  Earth  and  third-body  perturbations.  Consider  third-body 
perturbations  due  to  an  object  of  mass  M  at  a  geocentric  distance  D.  Since 
the  (ground  based)  radars  can  at  most  be  separated  by  2Re, 

c2AX/X  <  =  2GMRe/Dz 

-1 5 

When  the  third-body  is  the  Moon  AX/X  <  4.7  x  10  ,  when  it's  the  Sun 

AX/X  <  8.4  x  10'13. 

The  change  in  the  gravitational  potential  due  to  the  effects  of  the 
Earth's  oblateness  can  be  much  larger.  Using  the  expression 


-GMr 


1  - 


J2P2(sin<}>) 


14 


tor  the  geopotential  (J2  =  0.00108,  <p  =  terrestrial  latitude),  it  follows  that 


c2AX/X  <  2GMe  J2/Re3 
-9  . 

Hence,  AX/X  <  3.3  x  10  which  is  significantly  larger  than  the  second  order 
special  relativisitic  effect  due  to  the  Earth's  rotation  (oje  =  1 5V041  /sec , 

Re  =  6378.1  km  so  8E  =  0.47  km/sec/c  *  1.6  x  10~6,  62  =  2.4  x  10~12). 

The  second  general  relativistic  effect  is  due  to  the  lack  of  flatness 
(in  Minkowski  space)  of  the  general  relativistic  line  element  near  massive 
bodies.  This  results  in  a  departure  of  the  speed  of  light  near  such  an 
object  from  its  speed  in  vacuo  (eg.  c).  Hence  distance  to  an  object,  as 
deduced  by  halving  the  round  trip  flight  time  and  then  dividing  by  c,  must  be 
corrected.  In  the  lowest  order  post-Newtonian  approximation  the  amount  of 
the  correction  is  equal  to  the  path  integral  of  the  (Newtonian)  gravitational 
potential.  Thus, 

QU 

As  -  (2/c)  / Mt  «  (2/C)  /  -p  ^ 

“  (2GM£/c2 )  *n  (Rsat/RE) 


This  amounts  to  ^  1.7mm  for  a  near-stationary  satellite. 

There  are  still  two  approximations  within  the  context  of  a  special 
relativistic  computation  that  await  numerical  investigation.  One  is  =  6y 
and  the  other  is  =  -ftp  n^  =  -n^..  Consider  the  velocity  approximation 
first  when  the  transmitting  and  receiving  radars  are  colocated.  Then  the 
difference  between  8^  and  is  the  result  of  a  change  of  direction  due  to 
the  rotation  of  the  Earth  during  the  pulse's  travel  time.  Clearly  this  is 
of  order  A 8  =  uvAt  with  At  -  0?1  whence  we  are  cpeaking  of  an  effect  of  order 


10-5  (of 1  In  radians  Is  ^  10"^).  In  turn  the  lowest  order  (in  0E)  that 

this  correction  can  occur  is  first.  Hence  the  =  By  approximation  Is  of 
order  0£A6  ^  10“^8|.  Therefore,  in  any  analysis  good  to  62E,  this  term  can 
be  neglected.  An  analytical  refinement  of  this  argument,  directly  from  Eq.(7) 
or  Eq.  (11),  leads  to  the  same  result. 

The  other  approximation  is  the  neglect  of  an  aberration  difference,  n^ 
is  not  anti-parallel  to  because  they're  measured  in  two  different  inertial 

frames  (not  because  of  the  parallactic  displacement  of  the  radars).  But  the 
size  of  this  difference,  cf  Eq.  (1)  or  (2),  Is  of  the  order  of  the  velocity 
difference  between  the  two  reference  frames.  This  difference,  when  the  trans¬ 
mitting  radar  is  colocated  with  the  receiving  radar,  has  just  been  shown  to  be 
-5 

of  order  10  BE  in  magnitude.  As  the  unit  vectors  n^  and  ny  only  appear  in 
scalar  products  with  quantities  (at  most)  of  order  (10)  BE,  it  therefore 
follows  that  the  neglect  of  this  aberration  term  Is  of  order  10"5(10'4)B2e. 
Thus,  it  too  is  negligible  in  a  B2E  analysis. 

A  final  caution.  While  I  have  not  hesitated  to  evaluate  various  post- 
Newtonian  correction  terms  in  a  Newtonian  manner,  this  is  very  different 
from  looking  at  Eq.  (7)  and  deciding  that  an  ab  initio  Newtonian  derivation 
of  it  Is  appropriate.  That  would  be  wrong  and  can  clearly  never  lead  to  a 
post-Newtonian  term  to  evaluate.  Much  muddled  thinking  can  be  avoided  if  the 
problem  is  rigorously  solved  first  and  then  familiar  concepts  used  to  describe 
the  result. 
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